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1 Introduction 



Electromagnetic radiation from a uniformly accelerated charge is an example of an "eternal 
problem" of physics. Starting with the work by Born |]J different aspects of this problem 
have been discussed again and again (see e.g |2|, [3] and references therin). In this paper 
we consider a peculiar quantum analogue of this problem. 

Consider a small uncharged body and assume that it has internal degrees of freedom 
interacting with the electromagnetic field. A polarizable body is a well known example 
of such a system. Electromagnetic zero-point fluctuations induce dipole moments in the 
body. If a polarizable body is at rest, corrections to the field connected with this effect 
after averaging result in the change of the energy of the system. In the presence of two 
polarizable bodies the energy shift depends on the distance between them and results in 
the well known Casimir effect (see e.g. 0-0). For an accelerated motion of a polarizable 
body the net result of the emission of its induced fluctuating dipole moment is quantum 
radiation, or so called dynamical Casimir effect [|ITJ-[[nj]. For a non-relativistic motion and 
relativistic motion in 2 dimensional spacetime this effect is studied quite well, especially 
in the special case when the polarizability is very high and a surface of the body can be 



approximated by a reflecting mirror-like boundary ||8|— p2|. Much less results have been 
obtained for the relativistic motion in physical 4 dimensional spacetime. Exception are 
cases of a uniformely accelerated plane mirror |23| and relativistic expanding spherical 
mirrors |25j-[^6j. See also |27J where quantum radiation from an accelerated spherical 



body with a mirror boundary was considered. In this paper we study the effect of quantum 
radiation generated by an accelerated motion of a small polarizable body. 

To simplify calculations, we assume that internal degrees of freedom of a body interact 
not with the electromagnetic field but with a scalar massless field. We assume that 
a "refractive index" n inside a small ball differs from its vacuum value, 1. To solve the 
problem, we assume also that n is close to 1 and use the perturbation expansion in a small 
parameter n — 1. We consider a simplest accelerated motion when the direction and the 
value of the acceleration (as measured in the comoving frame) are constant. We calculate 
a correction to the vacuum Hadamard function created by a uniformly accelerated motion 
of the body. The main result of the paper is the expression for the quantum average energy 
density flux at infinity for this problem. 

The paper is organized as follows. Section 2 discusses the set up of the problem. We 
formulate the equation for a massless scalar field in a presence of "refracting" matter 
and discuss the case when a refractive ball is uniformely accelerated. Section 3 contains 
the calculation of the perturbed Hadamard function in the presence of the accelerated 
refractive ball of small size. We also discuss symmetry relations for observables at J + 
connected with the boost invariance of the problem. The expressions for (0 2 ) Icn and 
(T^) Icn in the wave zone are obtained and discussed in Section 4. Section 5 contains a 
discussion of the obtained results. 
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2 Model 



2.1 A "refractive" body in a static spacetime 

Our purpose is to study quantum radiation from a an accelerated "refractive" body. We 
shall use a simplified model by assuming that a body interacts with a quantum massless 
scalar field <p. 

In order to describe the model let us consider first a static gravitational field described 
by metric 

dS 2 = -A 2 dt 2 + 7^ dx* dx j . (2.1) 
Here A and 7^ do not depend on t\ £ M = d t is a Killing vector, and 

A 2 = -^^. (2.2) 



It is well known (see [28]) that the Maxwell equations in the metric (2.1) are identical 
to the Maxwell equations in a media with e = /i = A. Thus A' 1 plays the role of the 
refraction index n = 1/ y/eJI. We use this observation to introduce an effective refraction 
index into the equations for the scalar field. 

Consider a new metric dS 2 related to (|2.1| ) as follows 

A 2 

dS 2 = dt 2 + n^ij dx l dx^ , (2.3) 

n 

where the effective refraction index n = n(x % ) does not depend on time t. We choose an 
equation for the scalar field (p in the form 

n n( p = 0, (2.4) 

where D n is the "box" -operator in the metric ( |2.3| ). Equation 

dS 2 n = 0, (2.5) 



gives characteristics for the equation ( f2.4| ). Consider a Killing observer moving with a 
four-velocity = £ M /A In the reference frame of this observer, dr = Adt is a proper 
time, dl = yj 7^ dx % dxi is a proper distance, and equation ( |2.5|) takes the form 

(2.6) 

ar n 

The characteristics of the scalar field equation (|2.4| ) are rays moving with a velocity c/n 
with respect to a static observer, and hence n really plays the role of the refraction index. 
The equation (|2.4Q can be identically rewritten in the form 



n<p = D(n)(p, (2.7) 

where 

n 2 - 1 

D{n)<p=—g-%<p. (2.8) 
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In the case when [n — 1| <C 1, the term D(n)tp can be consided as a perturbation^. It should 
be emphasized that the form of the operator D(n) implies that equations ( p.7|) -( p^8|) can 
be easily generalized to the case of a media with dispersion. For a monochromatic wave 
of frequency u, d 2 — > — u 2 , and in order to take into account the dispersion it is sufficient 
to put n = n{uj). 

In our set-up we assume that n = 1 everywhere outside some four- dimensional region 
T where a refractive body is located, and n = const 7^ 1 inside this region. We assume 
that a body is static and rigid. Let B be a three-dimensional volume occupied by the 
body at t — to and S = dB be a surface of the body. Denote by E a three-dimensional 
surface formed by Killing trajectories passing through S. The region T is defined as a 
four-dimensional region located inside S. Thus we have S = dT. 

The operator fl2.8|) in a spacetime with such a body is of the form 

n 2 - 1 

D(n) = -g-m%> (2-9) 
where $ is the Heavyside step function. 

2.2 Uniformly accelerated body 

Now we adopt the above scheme to the case of a uniformly accelerated body moving in a 
flat spacetime. Let (T, X, Y, Z) be standard Cartesian coordinates so that the metric is 

dS 2 = -dT 2 + dX 2 + dY 2 + dZ 2 . (2.10) 

Denote by 7 a world line of a uniformly moving observer. If X is the direction of motion 
and a is the acceleration, 7 is described by the equation 

X 2 -T 2 = l 2 = a - 2 . (2.11) 

In what follows the length parameter / will play an important role. For this reason it is 
convenient from the very beginning to introduce dimensional coordinates 

T X Y Z , rt , 

i = T- X =T »=T z = l- (2 ' 12) 

Denote 

x = (1 + f) cosh 77, t = (1 + f)sinh77. (2.13) 

1 There is an ambiguity in the form of the metric ( |2.3| ). One can multiply the metric by any function 
f(n) which for n — 1 takes the value 1. This operation would modify the form of the operator D(n). 
In particular, a term proportional VnV would be generated. For a wave of characteristic frequency u> 
it gives a contribution ~ u/An/b where b is a size of the body. It can be considered as a perturbation 
only if it is much smaller than the leading derivative terms of the unperturbed operator which are of the 
order lo 2 . For our problem u> ~ a where a is the acceleration of the body, and ab <C 1. In order to escape 
problems connected with the applicability of the perturbation approach and to simplify the calculations 
we choose the special case / = 1. 



4 



Then the metric ( [2.10| ) takes the Rindler form 

dS 2 = I 2 [-(1 + O 2 dif + di 2 + dy 2 + dz 2 ] . (2.14) 
This metric is valid in the wedge x > \t\. The equation (|2.11|) for 7 takes the form 

^ = 0, (2.15) 

while If] is a proper time along 7. 

Surface 77 = rjQ is a plane with a three-dimensional flat metric. It is convenient to 
introduce spherical coordinates 0, 0) related to (£, j/, 2) as 

(£,y,z)=/3n\ (2.16) 

where 

n* = (cos 8, sin cos 0, sin 6 sin 0) (2-17) 

is a unit vector directed from the origin £ = y = z = 0to the point (£, y, z). 

Consider a small uniformly accelerated body, i.e. with the size much smaller than I. 
Such a body is at rest in the reference frame (|2.14j ). In the general case the surface of the 
body is described by the equation 

P = Po(e,v). (2.18) 



Relation fl2.18|) is the equation which defines the surface S, while fl2.18|) together with 



77 = 770 determines a position of the surface of the body at time rjQ. Later we consider a 
special case when the body is a ball and its surface is a sphere of radius b. For this case 

/3 = y = const. (2.19) 

3 Quantum Radiation of an Accelerated Refractive 
Body 

3.1 Hadamard functions and stress-energy tensor 

To find quantum radiation of an uniformly accelerated body with |n — 1| 1, we consider 
the operator D(n) as a perturbation and use perturbation expansion to obtain an answer. 
In the absence of the body the field equation is 

□y> = 0. (3.1) 

Denote 

G£\x,x') = (O\<p(x)0(x')+<p(x')<p(x)\O), (3.2) 
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the Hadamard function for the standard Minkowski vacuum |0). One has 

Gi 1) (x,x') = ^^^y, (3.3) 

where 

s 2 ( x , x ') = -{t- t') 2 + (x- x'f + (y- y'f + (z- z'f , (3.4) 

(t,x,y,z) and (i/, x', y', z') being dimensionless Cartesian coordinates of points x and x'. 
Consider the inhomogeneous equation 

n<p = j. (3.5) 

Its solution is 

tp(x) = <p (x) -I 4 1 G r *(x, x') j(x') dV . (3.6) 
where the retarded Green function Gq* is 

G?{x, aO = *(* a (*> *(* " • ( 3 - 7 ) 
By considering the right-hand side of Q2.7|) as a perturbation and using ( |3.6|) one gets 

(x, x') = G$ (x, x') + G« (x, x') , (3.8) 

where 

G«(x,x') = -Z 4 Ud A x" [G I Q et (x,x")D"G i o\x,x") 

+ G£V, x") D" G { o\x, x")] } . (3.9) 

Notation D" indicates that the operator D acts on the argument x" . We also use 
notation G^l for G^ — Gq. That is this object that is required for the calculation 
of physically observable quantities obtained by subtracting contribution of zero-point 
fluctuations in an empty spacetime. The integration in (|3.9|) is performed over the interior 
of the world-tube T. 

By calculating G^(x,x') we can find (<£> 2 ) ron and (T^) 1 ™ which characterize change in 
the fluctuations and in the stress-energy tensor generated by the moving body: 

{^{x))- = \ limGW(x,x'), (3.10) 

Z x —*x 

(7V(x)) rcn = lim G£l(x,x') . (3.11) 

x'— >x 

Here 

IV = Q - i) V (mV, + 9& 9 a ^ a V pi (v (/i V„) + V^V^) , (3.12) 

and £ is a parameter of non-minimal coupling. For £ = one gets the canonical stress- 
energy tensor, while for £ — 1/6 one gets the "improved" one with vanishing trace. 
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3.2 ((p 2 (x)) ren and (T^(x)) ren on J 



+ 



We shall perform calculations assuming that a point of observation, x, is located very far 
away from the moving body, in the so-called radiation zone. Denote 



1 x 2 + y 2 + z 2 , u — t — r, (x,y,z) — rN l , 

N l = (cos 6, sin 6 cos $, sin 9 sin $) . (3. 13) 

The coordinate u is the dimensionless retarded time, and (r, 0, $) are spherical coor- 
dinates in the inertial reference frame. The wave-zone corresponds to taking the limit 
r — > oo with k and iV* fixed. 

For calculations it is convenient to make the point splitting in t-direction, i.e. to choose 
x and x' (the arguments of G^H) so that 



r = r', N* = N* . 



We also put 



h 



U = Uq 



u = u + 



2 ' 2 

For such a point splitting, G^H(x, x') becomes a function of r, N l , uq and h: 



Gill(x,x') 



point 
splittcd 



G{r,u ,N\h). 



(3.14) 
(3.15) 

(3.16) 



bmce GB is a symmetric function of its coordinates x and x', G is an even function of 
its argument h. 

To obtain (y? 2 ) 10 ", ( p,10|) , it is sufficient to put h — in ( |3.16| ). The calculations of the 
energy density flux are more involved. One obtains 

dE 



u . _ R^Tuu)- = r 2 (T nu )" cn = 4vrr 2 lim G(r, u , iV\ /i) 



where 



n„ 



|(i-4e)al,-iflg 



(3.17) 



(3.18) 



As we shall see, the leading term of G in the wave zone is proportional to r 2 , so that 

1 



G(r,u ,N\h) 



r 2 l 2 



g l (u 0j e) + h 2 g 2 (u ,e) + o(h 4 



(3.19) 



We include factor I 2 to restore the correct dimensionality of G. Functions Q\ and Qi are 
dimensionless. They do not depend on $ since the system is invariant with respect to 
rotation in (y,z)— plane. Combining these results, we get 

Gi(u , 9) 



2r 2 l 2 



dE 1 
dUdn I 2 



-(l-^^Kei-^Ke) 



(3.20) 
(3.21) 
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3.3 Boost invariance 



An important additional information about (ip 2 ) rcn and the energy flux at J + can be ob- 
tained by using the symmetry of the problem. The spacetime with a uniformly accelerated 
body is invariant under boost transformations 

t — > t = j(t + vx) , 

x^x = j(x + vt), (3.22) 

V = V, z = z } 7 = (1 -v 2 y l/2 . 

It is easy to determine how this symmetry transformation acts on J + . To do this, we 
introduce retarded spherical coordinates (u, r, 0, $) and (u, f , 0, $) for both Minkowski 
frames x and x, and using ( [3.22 ), we find the relation between them. In the wave-zone 
limit, r — > oo, (u, 0, $) - fixed, we have 

vu cos 



r = jr(l + v cos 0) + 7 



l + v cos 



u 



r 



tan 



l + v cos ' 
sin 



(3.23) 



$ = $ 



7(cos + v ) 

The invariance of (y? 2 ) 1 '™ under the transformation ( p. 22 ) requies that 

Gi(u,Q) =7 2 (1 + ^cos0) 2 6? 1 (m,0). 



(3.24) 



The invariance condition ( 3.24 ) can be written in the infinitesimal form. For this we put 
v = 5v and put the first variation of ( |3.24 ) with respect to 8v at v = equal to zero. This 
gives the following relation 



du 



'd®\ 



-7^) d u g 1 +\ — \ d e g 1 = 2coaeg 1 . 



Using ( p.23| ), we get 



dv 



du 
dv 



v=0 



dv I „ 

/ v=0 



(3.25) 



—u cos 



v=0 



'd& 

dv 



— sin . 



Hence we have 



dg 1 



ggi | 

d(\nu) (91n(sin0) 



v=0 



-2G 



1 ■ 



A general solution of the equation ( 3.27]) can be presented in the form 

Gi(u,e) 

In a similar way we get 



lysine 



u 



1 sin©, 
— ri2\ 



u 



(3.26) 
(3.27) 

(3.28) 
(3.29) 
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4 ((^ 2 ) ren and (T^) ren in the Wave Zone 
4.1 Wave zone approximation 

To calculate functions Q\ and Q2, we rewrite ( |3.9|) in a more explicit form. Since integration 
is performed over the world tube T, it is convenient to write the integral in ( |3.9| ) in the 
Rindler coordinates 

j^x" ...= J d V J du j^ {6,4> (3 2 Ad{3 ...= J d\ ... , (4.1) 

where A — 1 + j3 cos 9, du = sin 9 d9 d<p and (3q defines the boundary of the body, see 
( P-20| ). In these coordinates we also have 

(n 2 -1)V 

D = 1 p ) , (4.2) 

V = A- 2 d{fo{9^)-P) d 2 . (4.3) 

Both Green functions G T ^ and depend only on the distance s between a point 
in the wave zone and a point inside or on the boundary of the world tube T. Simple 
calculations give 

s 2 (x, x") = -2r w -u 2 + 2uA sinh rj + 1 + 2f3n l + (3 2 , (4.4) 

where 

w = u + (3 n±N± — F(rj) A , F(rj) = sinh 77 — iV 1 cosh 77 . (4.5) 

Here (u,r, 0, $) are retarded spherical coordinates of the point x in the wave zone, and 
(7], (3, 9, <fi) are Rindler spherical coordinates of the point x" in the tube T. The vectors 
n and AT" are defined by equations (|2.17|) and ( |3.13|) , respectively, and n± = (0,n 2 ,n 3 ), 
N ± = (0,N 2 ,N 3 ). 

An important observation is that the leading, 1/r 2 , term of G^H in the wave zone 
can be obtained by neglecting the terms independent of r in ( f4.4| ), that is by using the 
following approximate expressions 

s 2 (x,x") ~ — 2rw- , s 2 (x',x") ~ —2rw + , (4.6) 

where 

h 

w± = w ±-, w = u + pn ± N L - F{r,) A. (4.7) 
Combining all these results and using (|3.19|) , we get 

gi{u ,e) + h 2 g 2 { Uo ,e) + ... 



n 2 -l , 

167T" ■' lr 
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8{w_)V{— ) + 5(w+)T>(— ) 



(4.8) 



We omitted ^-function which enters the definition ( |3.7|) of Gq * since a future-directed null 
cone emitted from a point in the wave zone never crosses the tube T. 
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4.2 Calculation of integrals 

Denote by I(rjo, h) the following integral 

/OO 1 
d V 6(w-)%{—). 
-oo ' W + 

Notice that 



FA 2(F') 2 A 2 



1 w+ 



w± 



+ 



w± 



where ( )' = d„( ). We use here the following property of the function F(rj) 



(4.9) 
(4.10) 

(4.11) 



For 0/0, F(r]) is a monotonically increasing function which changes from — oo (at 
rj = — oo) to +oo (at rj = oo). Thus equation w (r]) = c has a unique solution for any c. 
We denote by rj and r)± the solutions of the following equations 



The integral ( |4.9| ) can be easily calculated by using the relation 

Kv-) 



F ; (r)-)A 



We get 



F'{r}- 



Now we use the following Taylor expansions 



oo fan 

v- = vo + 22 t Zn ' 



w (V-) 



n=l 

oo 



111 



AT,^Fi n) (V--VoT 



n=l 



111 



(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 



where F {n) = F^(r) ). Substituting Q435D into fl4T| ) and using (|4.11|) , we solve equation 

MV-) = \ (4-17) 
to determine z n in terms of Fq and Fq. For calculations we use Maple. The result is 

1 F _ 3 (F ) 2 + (^) 2 



2AFL ' 



22 



2 4 



4y4 2 (F(5) 3 ' 



^3 



8^ (F ') 



3F [-5(Fq) 2 + 3( ^ 
16A 4 (I^) 7 



A2l 



(4.18) 



z 5 



3[35(F ) 4 -30(F ) 2 (^) 2 + 3(^) 4 ] 



32A 5 (i^) 9 

We performed calculations up to the order h 5 which is required to calculate fl4.9|) up to 
the order h 2 . 

Next steps of the calculations are the following: 
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1. Write F(?7_) as 

F(r]J) = F cosh(r]_ - ?7o) + Fg sinh(r]_ - r] ) 



(4.19) 



2. Substitute ( |4.15| ) into expression ( |4.19| ) and using ( [4. 18| ) get the expression for F(rj_) 
in powers of /i. 



3. Use the obtained expansion to calculate I(r]o, h) defined by ( |4.14j) . 

4. Calculate Jo(r] , h) = I(r] , h) + Jo^o, — h) to obtain the expression for 



J{Vo,h) 



ch] 



5(w„)d 2 (—) + 5(w + )d 2 / ' 



(4.20) 



We performed these calculations using Maple. The final result is 



J{Vo, h ) 



F sin 2 6 



1 + 



h 2 



8A 2 (F^ 



— [(4 cosh 2 7] + 3) cos 2 



2A*(F>y 

— 8 sinh r] cosh r] cos — 7 + 4 cosh 2 r/o] 
Thus the right-hand side of ([4.8|) is 

J du ^ dp^A- 1 J( Vo ,h). 
In this relation r] is a solution of the equation 

F(r ]0 ) = A- 1 (u + (3n ± N ± ). 



(4.21) 



(4.22) 



(4.23) 



Since < 1, one can solve this equation perturbatively. Let 770 be a solution of the 
equation 

F(f, ) = u , (4.24) 



then 



770 - ^0 ~ r=r- {n ± Nj 
r 



U COS I 



(4.25) 



where F' = F'(fj ). In order to take into account the dependence of J(r) ,h) on (3 it 
is sufficient to expand J near fjo and use ( 4.25|) . The obtained corrections contain an 
additional small factor b/l, where b is the size of the body. By neglecting this correction 
and similar corrections in A, we get 



in 2 -1)V 



Jo{fjo,h), 



16tt 3 P 

where V is the volume of the body. By comparison with Q4.8| ) we get 

F sin 2 



g 1 (u ,e) = -B- 



in) 



(4.26) 



(4.27) 
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^( M o,e) = -S^^(4F 2 -3sin 2 (e)) , (4.28) 

where 

« (n 2 -l)V , , 

B=L l2^- (429) 
Let us emphasize once again that since we are considering the leading in b/l terms, to 

calculate F and Fq one must put 6 = 0. In particular, F = u . In order to get Fq one 
must first solve the equation 

F = sinh ?7o — sin 6 cosh r] = u , (4.30) 

and determine r] Q = T]o(uo, cos0), ans substitute this value into the definition of Fq 

Fq = cosh r] — sin 6 sinh r] . (4-31) 



4.3 (cp 2 ) ren and energy density flux 



It is easy to check that Qy^ obey the symmetry relations ( |3.28| )- fl3.29| ). Really by using 
relations ( |4.30| ) and ( [4.31| ) one can rewrite in the form 



Q x [u,Q) = -B^-, g 2 {u,e) = -B^, (4.32) 



where z = sin(0)/w and 



z 2 z 2 U-3z 2 ) 

^= (1 + 3 2)5/2 ' M = 8( 1 + ,2)9/2 ■ ( 433 ) 



Using equations ( |3.20| ) and ( |3.21| ), we get 



. 0l (n 2 — l)6 3 aqi (z) 

(iE 1 (n 2 — l)b 3 a . , , , 

^1 -50^)- (4-35) 



Here 

1 - 3z 2 /4 

In the previous relations we restored dimensional coordinates R and U = T — R, and z 
which enters these relations is sin(0)/(ai7), where a is the acceleration of the body. Plots 
of functions gi and #3 are shown in Figure |T] 

By integrating over angles one gets the total energy density flux 

where 

, 1 + 4u 2 + 5u 4 + 10m 6 

/(») = • ( 4 - 38 ) 
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Figure 1: Plots of functions gi(z) (gl) and gs(z) (#3). 



5 Conclusion 

In this paper we considered the vacuum polarization effects in the presence of a uniformly 
accelerated dielectic body. We considered a scalar field model. Under assumption that 
a size of the body, b, is much smaller than the inverse acceleration, a -1 , we calculated 
the energy density flux created by the accelerated body at infinity. This flux is given by 
( 4.37[ ), and for the canonical energy (i.e. for £ = 0) it is always negative. Its divergence 



~ U~ 2 near U — is connected with the idealization of the problem: it is assumed 
that the motion remains uniformly accelerated for an infinite interval of time. The boost- 
invariance property connected with this assumption significantly simplifies calculations. In 
particular, the quantity dE / (dU dVt) giving the anglular distribution of the energy density 
flux at J + at given moment of the retarted time U, besides common scale dependence, 
£7~ 4 , depends on the function of one variable, sin Q/ all. 

It would be interesting to repeat calculations for a more realistic case of the electromag- 
netic field. One can expect that some general features, especially those that are related to 
the symmetry of the problem, will remain similar to the case of the scalar massless field, 
while details, such as the angular distribution of the energy flux which might depend on 
the spin of the field will differ. It should be specially emphasized that the dependence of 
the electromagnetic field equations on the dielectric and magnetic properties of the media 
is uniquely fixed, while in the model case of a scalar field there is an ambiguity which we 
specially fixed to simplify the model. 

In our consideration we assumed that a uniformly accelerated refractive body is cold, 
that is its temperature is zero. One can also consider a case when a body is heated. 
Especialy interesting is a case when the temperature of he body coincides with the Unruh 



temperature corresponding to its acceleration. General arguments given in |29], |30| allows 
one to expect that in this case the quantum radiation vanishes. We hope to return to 
these problems somewhere else. 
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